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ABSTRACT 

The thermal and fragmentation properties of star-forming clouds have important con- 
sequences on the corresponding characteristic stellar mass. The initial composition of 
the gas within these clouds is a record of the nuclcosynthctic products of previous 
stellar generations. In this paper we present a model for the evolution of star- forming 
clouds enriched by metals and dust from the first supernovae, resulting from the explo- 
sions of metal- free progenitors with masses in the range 12 - 30 Mq and 140 - 260 Mq. 
Using a self-consistent approach, we show that: (i) metals depleted onto dust grains 
play a fundamental role, enabling fragmentation to solar or sub-solar mass scales al- 
ready at metallicities Zcr — IQ^^Zq] (ii) even at metalhcities as high as 1Q~'^Zq, 
metals diffused in the gas-phase lead to fragment mass scales which are ^ IOOM0; 
(iii) C atoms are strongly depleted onto amorphous carbon grains and CO molecules 
so that CII plays a minor role in gas cooling, leaving 01 as the main gas-phase cooling 
agent in low-metallicity clouds. These conclusions hold independently of the assumed 
supernova progenitors and suggest that the onset of low-mass star formation is con- 
ditioned to the presence of dust in the parent clouds. 

Key words: stars: formation, population III, supernovae: general - cosmology: theory 
- galaxies: evolution, stellar content - ISM: abundances, dust 



1 INTRODUCTION 

Supernovae represent the first cosmic polluters. On 
timescales shorter than 50 Myr, these stellar explosions seed 
the cosmic gas with the first metals. Theoretical models sug- 
gest that a significant fraction of these heavy elements could 
be in the form of dust grains (Kozasa, Hasegawa & Nomoto 
1989; Todini & Ferrara 2001; Nozawa et al. 2003; Schnei- 
der, Ferrara & Salvaterra 2004). Indeed, following the ex- 
plosion, conditions are met in the expanding SN ejecta for 
the condensation of gas-phase metals in grains of different 
size and species. In particular, the fraction of metals de- 
pleted onto dust, /dcp = Afdust /A/mct , varies between 0.3 
and 0.7 for pair-instability supernovae (PISN) with progeni- 
tor masses in the range 140 — 260Mq (Schneider et al. 2004), 
and between 0.24 and 1 for Type-II supernovae (SNII) with 
12 - 3OM0 progenitors (Todini & Ferrara 2001). 

Efficient dust production on short timescales in su- 
pernova ejecta provides a straightforward interpretation of 
the large dust masses (Mdust '~ 10^ Mq) inferred from mm 
and submm observations of quasars (QSOs) with redshifts 
5.7 < z < 6.42, (Bertoldi et al. 2003: Priddey et al. 2003). 



Dust condensation in the atmospheres of AGB stars requires 
timescales of 1 - 2 Gyr (Morgan & Edmunds 2003), too long 
to account for the presence of large dust masses at redshifts 
z > 5, when the age of the Universe is < 1.2 Gyr. Fur- 
ther support to a SN origin for dust at high redshift has 
recently come from near-IR observations of a z = 6.22 QSO 
(Maiolino et al. 2004). The extinction properties of the ob- 
served spectrum cannot be reproduced assuming a Small 
Magellanic Cloud dust extinction curve (which generally ap- 
plies to z < 5 systems) and can, instead, be interpreted with 
models for dust produced in SNII ejecta. Still, direct IR ob- 
servations of SN remnants (Cassiopea A, Crab, Kepler) have 
detected only thermal emission from a hot dust component 
(temperatures 50 — 80 K) with associated dust masses of 
only 0.003 - O.O7M0 (Hines et al. 2004; Green et al. 2004). 
A larger dust component could be at lower temperatures 
(< 30 K), and thus not visible at IR wavelengths. Submm 
observations of some of these systems have, so far, yielded 
only upper limits, with estimated masses < 0.2 — I.^Mq 
(Morgan et al. 2003; Krause et al. 2004; Wilson & Batrla 
2005). 

Early metal enrichment of the cosmic gas has important 
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implications for the formation of the first low-mass and long- 
lived stars. In the last few years, results of sophisticated nu- 
merical simulations (Abel, Brian & Norman 2002; Bromm, 
Coppi & Larson 2002) and semi-analytical models (Omukai 
& Nishi 1998; Omukai 2000; Nakamura & Umomura 2002) 
have shown that, due to the primordial composition of the 
gas, the first stars were predominantly very massive, with 
masses of 100 — 1000 Mq. The presence of metals favors 
cooling and fragmentation of prestellar-gas clouds, enabling 
the formation of solar and sub-solar mass stars when the gas 
initial metallicity reaches a critical value Z^t- Despite its im- 
portance, theoretical studies have only constrained this fun- 
damental parameter to within a few orders of magnitudes, 
W^^Zq < Zci < 5 X 10"'^ Z0, with the lowest limit implied 
by models which include the effect of dust cooling (Schnei- 
der, Ferraxa, Natarajan & Omukai 2002; Schneider et al. 
2003; Omukai, Tsuribe, Schneider & Ferrara 2005) and the 
highest limit suggested by studies where metal line cooling is 
considered to be the main driver of the transition (Bromm, 
Ferrara, Coppi & Larson 1999). 

All these theoretical studies relics on the assumption 
that the relative metal abundances and dust grains proper- 
ties which characterize the initial composition of pre-stellar 
gas clouds are the same as those observed in the local in- 
terstellar medium, simply rescaled to lower values of metal- 
licity. However, the composition of the cosmic gas enriched 
by nucleosyrithetic products of the first stellar generation 
through supernova explosions is likely to be different. To 
overcome this basic limitation, Bromm & Loeb (2003) have 
derived individual critical abundances of C and O but have 
limited their analysis to CII and 01 fine-structure line cool- 
ing during the relatively low-density evolution of collapsing 
pre-stellar clouds. 

In this paper, we follow a different approach and study 
the detailed thermal and fragmentation history of collapsing 
clouds enriched by the first supernova explosions. We use the 
model developed by Omukai (2000), and recently improved 
in Omukai et al. (2005) to carry out the first self-consistent 
study which explicitly accounts for the relative metal abun- 
dances and dust properties predicted by supernova models. 
We follow the fate of collapsing clouds enriched by these nu- 
cleosynthetic products to investigate whether preferred mass 
scales can be identified, and eventually linked to character- 
istic masses of second generation stars. 

The initial composition of collapsing gas clouds must be 
inferred by results of stellar and supernova models of metal- 
free stars. This relies on specific assumptions on the initial 
mass function (IMF) of Population III stars. Unfortunately, 
both the mass range and the shape of the IMF are still 
hampered from our limited understanding of the accretion 
physics and protostcUar feedback effects (Omukai & Palla 
2003; Bromm & Loeb 2004; Tan & McKoo 2004; see also 
Bromm & Laxson 2004 for a thorough review) . Given these 
uncertainties, we consider Population III stars forming both 
in the SNII and PISN progenitor mass ranges, taking 22 Mr.-, 
and 195 Mq stars as representative cases for the two classes. 

The paper is organized as follows: in section 2 we de- 
scribe the characteristics of the model, with particular at- 
tention on those aspects which have been specifically intro- 
duced to account for the properties of supernova dust grains. 
In section 3 we present the results of the analysis, and in sec- 



tion 4 we summarize the main conclusions and discuss their 
implications. 



2 THE MODEL 

The thermal and chemical evolution of star forming gas 
clouds is studied using the model developed by Omukai 
(2000) and recently improved in Omukai et al. (2005). The 
collapsing gas clouds are described by a one-zone approach 
where all physical quantities are evaluated at the center as a 
function of the central density of hydrogen nuclei, uh- The 
temperature evolution is computed by solving the energy 
equation 

de d 1 . ... 

dt at p 

where the pressure, p, and the specific thermal energy, e, are 
given by, 

pkT 



P = 



prriH 
1 kT 



(2) 



(3) 



7 - 1 /i TTlH 

and p is the central density, T is the temperature, 7 is the 
adiabatic exponent, p, is the mean molecular weight and mn 
is the mass of hydrogen nuclei. The terms on the right-hand 
side of the energy equation are the compressional heating 

rate, 



tg = 



31T 

32Gp' 



(4) 



and the net cooling rate, Anet , which consists of three com- 
ponents. 



Anct — Alino + Acont + Ac 



(5) 



The first component, Anne, represents the cooling rate due 
to the emission of line radiation, which includes molecular 
line emission of H2, HD, OH and H2O, and atomic fine- 
structure line emission of CI, CII and 01. Following Omukai 
et al. (2005), H2 coUisional transition rates and HD parame- 
ters are taken from Galli & Palla (1998). The second compo- 
nent, Acont, represents the cooling rate due to the emission 
of continuum radiation, which includes thermal emission by 
dust grains and H2 collision-induced emission. The last term, 
Achom, indicates the cooling/heating rate due to chemical 
reactions. When the gas cloud is optically thick to a spe- 
cific cooling radiation, the cooling rate is correspondingly 
reduced by multiplying by the photon escape probability. 
Unless otherwise stated, the treatment of these processes 
are the same as in Omukai (2000), to which we refer the 
interested reader for further details. 

In the present analysis, we relax the assumptions that 
gas-phase elemental abundances and dust grain properties 
are the same as those observed in local interstellar clouds. 
Instead, we compute the relative abundances of gas-phase 
metals and dust grains applying the model of Todini & Fer- 
rara (2001) for metal-free Type-II supernovae (SNII) and of 
Schneider, Ferrara & Salvatcrra (2003) for pair-instability 
supernovae (PISN). The properties of dust (relative dust 
compounds and grain size distributions) affect two processes 
which contribute to the thermal evolution of gas clouds: (i) 
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the cooling rate due to dust thermal emission and (ii) the 
H2 formation rate on dust grains. The description of these 
processes has been modified with respect to the treatment 
in Omukai (2000) and Omukai et al. (2005) to properly take 
into account the properties of dust grains produced in the 
first supernova explosions. 



2.1 Dust cooling 

The contribution of dust grains to gas cooling is determined 

by the balance between dust thermal emission and heating 
due to collisions with gas particles, 



4 a Tgr Kp Pe: 



where a is the Stefan-Boltzmann constant, Tg 



(6) 

is the dust 



temperature, Kp is the Planck mean opacity of dust grains 
per unit dust mass, /Jesc is the photon escape probability, 
and the dust mass density pgr can be expressed as pV, with 
D indicating the dust-to-gas ratio. Following Omukai (2000) 
the photon escape probability can be written as. 



:min(l,T ) 



with 



KppXj 



(7) 



where the optical depth is estimated across one local Jeans 
length (roughly the size of the central core region in the 
Penston-Larson similarity solution). The Planck mean opac- 
ity per unit dust mass is, 



Kp 



B„{Tgr)K{u)dl' 



(8) 



where k(i') is the absorption coefficient for the frequency u 
per unit dust mass and By(Tgr) is the Planck brightness. 

The right-hand side of the eq. (6) is the coUisional heat- 
ing (cooling) rate per unit volume (erg s~^ cm~^) of dust 
(gas) particles and can be written as (HoUenbach & McKee 
1979), 



tcoll 



(9) 



where = nu Cgr vh f is the average time between two 
successive collisions, rigr and dgr arc the grain number den- 
sity and cross section, vh is the average speed of hydrogen 
nuclei and / — nv/nHVH measures the contribution of other 
species. We assume that the gas is fully molecular and fol- 
lows a Maxwellian distribution so that. 



= (- 
Vtt? 



, , (10) 

and / = 0.3536 -|- 0.5 yue, where yue = riHe/nH = 0.083 
(which corresponds to a He mass fraction of Y=0.25) and we 
have neglected the contribution of heavy elements. These are 
good approximations in the density range where collisional 
dust heating is relevant (nn > 10^^~^^cm~^) and for the gas 
metallicities that we are interested to {Z = Za). 
If we further write, 



nn mH(l + 4j/ho) SV 



(11) 



in eq. (9), where S is the total grain cross section per unit 
mass of dust, eq. (6) appears to be independent of the dust- 
to-gas ratio and to depend on the thermal states of gas 
and dust and on grain specific properties (cross section and 
Planck mean opacity). 
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Figure 1. Absorption coefficients per unit dust mass assuming 
the properties (size distributions and relative masses of different 
grain species) of dust grains predicted to form in the ejecta of a 
195 Mq PISN (solid lines) and of a 22 M© SN (dashed lines). 



Given the grain size distribution functions per unit dust 
mass (cm^"'^ gr^"*^), ngr_i(a) = c?ngr,i(a)/(io, and the mass 
fractional abundance, of each grain type i predicted by 
the SN-dust formation models, the total grain cross section 
per unit dust mass (cm^/gr) is. 



5 = Ei /i Si with Si 



Jo 



ngr,i(a)7ra da. (12) 



Similarly, the total absorption coefficient per unit dust mass 

(cm^/gr) is. 



k{i') = Ei fi with 



/•oc 

Kii^i — I Q 1,(^(1 

Jo 



) ngr,i(a) TTO rfo,(13) 



where Ql are the absorption (or extinction) cross section 
normalized to the geometrical cross-sections. For each dust 
species, the Q parameters for any frequency and radius have 
been computed using the standard Mie theory for spherical 
grains with grain optical constants collected from published 
data (see Table 1 for details). In Fig. 1 we show the to- 
tal absorption coefficients per unit dust mass assuming the 
properties of dust grains formed in the ejecta of a 195 Mq 
PISN (see Figs. 4 and 5 in Schneider, Ferrara & Salvaterra 
2004) and of a Z=0, 22 Mq SN (see Figs. 5 and 6 in To- 
dini & Ferrara 2001). For the PISN model, the absorption 
coefficient is dominated by SiOg grains in the 6 /im - 200 
fiui range, whereas above 200 /im amorphous carbon (AC) 
grains give a comparable contribution. Sliortward of 6 pm, 
AC grains largely dominate down to ~ 0.2/im, below which 
Mg2Si04 grains produce the steepening of the curve. For the 
SNII model, the absorption coefficient is largely dominated 
by AC grains except for the 0.8 pm - 30 pm range, where 
the Mg2Si04 features appear. 

The resulting Planck mean opacities are shown in Fig. 2. 
The solid line shows the absorption Planck mean per unit 
dust mass for the 195 A4q PISN model and the dashed line 
shows the same result for the 22 Mq model. When com- 
puting eq. (8), it is important to take into account the ef- 
fect of dust sublimation: grains can be destroyed by heating 
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Table 1. Properties of dust grains produced in the ejecta of a Z=0, 22 Mq SNII and of a 195 Mq PISN. The table show various grain 
species (first column), their mass fraction relative to the total dust mass produced (second and fourth columns), and their geometrical 
cross-sections (third and fifth columns). References for the adopted grain optical properties are given in the last column. In some cases, 
data was linearly extrapolated on the log-log plot to extend the original wavelength coverage to longer wavelengths. 
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Figure 2. Planck mean opacities per unit dust mass for absorp- 
tion for the 195 Mq PISN model (solid) and for the 22 Mq SN 
case (dashed). Thin lines assume that the grains never sublimate 
(see text). 



through collisions with gas particles if the equilibrium grain 
temperature Tgr exceeds a characteristic value Tl^^ which, 
for each grain species i, depends on the gas density through 
the relative partial pressures of the heavy elements in the 
gas phase. Here we assume a sublimation temperature of 
= 2000 K for AC grains and of 1500 K for other grain 
species (Laor & Draine 1993; Preibisch et al. 1993) and ig- 
nore the density dependence *. The thin lines in Fig. 2 show 
the results when sublimation is not considered. 



electron shells of H and surface atoms (van der Waals in- 
teractions), and a chemisorbed site, when overlap between 
their respective wave-functions occur. Coming from the gas 
phase, the H atoms are first physisorbed and then either 
cross the barrier between the physisorbed and chemisorbed 
sites (by moving perpendicularly to the surface) , or go to an- 
other physisorbed site (by moving along the surface). The 
relevant parameters which regulate grains surface character- 
istics are therefore the desorption energies of physisorbed H, 
Eup, of chemisorbed H, Euc, and of the saddle point be- 
tween a physisorbed and a chemisorbed site, Es, and can be 
derived by laboratory experiments for different grain species 
(see Cazeaux & Spaans 2004 and references therein) 1^. 

At temperatures where the process of H2 formation on 
grain surface is relevant in collapsing gas clouds, the most 
efficient reaction is the "collision" of two chemisorbed atoms, 
with efficiency (Cazeaux & Tielens 2002), 



<:h2 = 1 + 



QPC 



(14) 



where we have neglected the correction factor at high tem- 
peratures (T > 300 K) which reflects the evaporation of 
chemisorbed H. In the above equation. Pup is the evapora- 
tion rate of physisorbed H and apc is the mobility to go 
from a physisorbed to a chemisorbed site, and their ratio 
can be approximated as (Cazeaux & Tielens 2002), 



(15) 




V 


Es " 


exp 





The resulting H2 formation rate per unit volume (cm 
on grain surface can be expressed as. 



2.2 H2 formation on dust grains 

An H atom can be bound to the grain surface in two energet- 
ically different sites: a physisorbed site, when the attractive 
forces are due to mutually induced dipole moments in the 



-Rh2 = 2 '^g'^ '^g'^ '=^2 'S'h (T) 

= fcgr n(H) nu 



(16) 



where n(H) is the number density of H atoms in the gas 
phase, riH is the number density of H nuclei (protons), and 



* Note that comparable values are found for the condensation 
temperatures in the 195 Mq model, namely 2100 K for AC grains, 
1800 K for AI2O3, and 1470 K < T < 1560 K for the other grain 
species (see Todini & Ferrara 2001; Schneider et al. 2004). 



t From Table 1 of Cazaux & Spaans (2004), we take i?HP = 
800 K, Ehc = 3 X 10"' K, and Es = 250 K for AC grains. For other 
grain species we assume the same parameters as those derived for 
olivine (Mg2Si04), namely Eup = 650 K, Euc = 3 X 10"" K, and 
Es = 200 K 
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Sh{T) is the sticking coefiicient of H atoms which depends 
both on dust and gas temperatures (HoUenbach & McKee 
1979), 

5h(T) = (17) 
1 + 0-4 (^) ' + 0.2 (j^) + 0.08 {j^y . 

Substituting eqs. (10) and (11) in eq. (16), the coefficient 

fcgr can be written as, 

fegr = f^)'^' SH{T)VmH{l + 4yue)SifiSieii^, (18) 

where is the H2 recombination efficiency for each grain 
species. 



3 RESULTS 

3.1 Thermal evolution of collapsing clouds 

Depending on the efficiency of metal diffusion and mixing, 
gas clouds in the vicinity of supernovae can be enriched to a 
wide range of metallicity values. If we assume that gas-phase 
metals and dust grains are not selectively transported, we 
can take the initial composition of pro-stoUar gas clouds to 
reflect the final composition of SN ejecta. We then follow 
the thermal evolution of collapsing gas for a set of initial 
metallicitios but relative metal and dust properties consis- 
tent with those predicted by the corresponding supernova 
model. 

The results are shown in Fig. 3 for gas clouds with ini- 
tial metallicities in the range < Z/Zq < 0.1 enriched 

by the products of a 195 Mq PISN (left panel) and of a 
22 Mq metal- free SNll (right panel). For these specific su- 
pernova models, dust depletion factors are predicted to be 
/dep = 0.65 and 0.24, respectively. Consistent with previous 
findings (Omukai 2000; Schneider et al. 2002, 2003; Omukai 
et al. 2005), at relatively low densities (nn ^ lO^cm"^), 
thermal evolution is dominated by line-cooling of H2 for 
Z/Zq < 10-'', of HD for 10"* < Z/Zq < 10-^ and, for 
Z/Zq ^ 10^^, of O fine structure lines. Once the relevant 
cooling agent reaches the NLTE - LTE transition or the gas 
becomes optically thick to the dominant cooling radiation, 
line-cooling becomes inefficient and gravitational contrac- 
tion leads to a temperature increase. 

Note that, contrary to what commonly assumed, the 
above results suggest that C makes a negligible contribution 
to fine-structure line cooling, being strongly depleted onto 
AC grains and CO molecules for most SN models (Todini & 
Ferrara 2001; Schneider et al. 2004). 

The impact of dust-cooling on the thermal evolution 
starts to be apparent for initial metallicities Z/Zq ^ lO^'' 
{Z > Zci) as a temperature dip at high densities 
(riH ^ lO'^^cm"^), which progressively shifts to lower densi- 
ties for increasing total metallicity^ , until, for Z/Zq ^ 10~^, 

■t If a larger amount of dust is present because the total metallic- 
ity of the gas is higher, then the corresponding cooling rate starts 
to be effective earlier in the evolution of collapsing gas clouds and 
it is not confined to the highest density regime. 



the temperature dips due to line-cooling and dust-cooling 
merge and it is no-longer possible to separate the two con- 
tributions. 



3.2 The impact of SN dust 

It is interesting to compare the results of the self-consistent 
model for thermal evolution of collapsing clouds with pre- 
vious studies, such as in Omukai et al. (2005), based on 
different assumptions on dust and metal properties. Com- 
paring the tracks shown in Fig. 3 with the results in Fig. 1 
of Omukai et al. (2005), we see that the temperature dip 
at high density (nu ~ 10^^ cm~^) due to dust-cooling 
which in the present analysis is already apparent in the 
Z = W~^Zq track, starts to be evident only at higher 
metallicities {Z > lO'^'Ze) in Omukai ct al. (2005). This 
result reflects a higher cooling efficiency of SN-dust with re- 
spect to present-day interstellar dust, partly due to larger 
grain cross-section and absorption coefficient per unit dust 
mass (eqs. 12-13). Indeed, due the limited time available for 
grain-growth, SN-dust grains are smaller than present-day 
interstellar grains and the ratio of area to mass tends to be 
larger for smaller grains (Todini & Ferrara 2001; Nozawa et 
al. 2003; Schneider et al. 2004). A major difference is repre- 
sented by the condensed C species adopted: in SN-dust mod- 
els, C atoms condense in AC grains which arc characterized 
by a high sublimation temperature (see section 2.1) and sur- 
vive to the high-temperature environment of low-metallicity 
clouds. Conversely, following Pollack ct al. (1994) all pre- 
vious models assumed that C atoms condense in organics 
and therefore sublimate at lower temperatures, ^ 600 K 
(Omukai 2000; Schneider et al. 2003; Omukai et al. 2005). 

The larger grain cross-section per unit mass of SN-dust 
enhances the H2 formation rate. However, this has only a 
minor effect on the thermal evolutionary tracks at the lowest 
metallicities, being relevant only for Z ^ W~'^Zq. 



3.3 Fragmentation properties 

Our main interest here is to investigate the fragmentation 
properties of clouds enriched by primordial supernovae, with 
the aim of understanding whether a preferred mass scale 
can be identified and eventually related to the characteris- 
tic mass of second-generation stars. Indeed, numerical sim- 
ulations of primordial star formation have shown that in 
primordial environments, where magnetic fields, turbulence 
and rotation are likely to play a minor role in the dynamics, 
a preferred mass scale can be related to the Jeans mass at 
the end of the fragmentation process (sec Bromm & Larson 
2004 and references therein). More generally, Larson (2005) 
suggests that the low-mass, Salpeter-like IMF, which char- 
acterizes contemporary star formation, is determined largely 
by thermal physics and fragmentation in pro-stellar clouds. 

Fragmentation occurs efficiently when the adiabatic in- 
dex 7 — dlogp/dlogp < 1, and almost stops when isother- 
mality breaks (7 > 1) as also shown by the simulations of Li, 
Klesser & MacLow (2003). Thus, consistent with Schneider 
et al. (2002), (2003) wc can adopt the condition 7 = 1 to 
identify the preferred mass scale of the fragments (Inutsuka 
& Miyama 1997; Jappsen et al. 2005) 

In Fig. 4 we plot the pre-stellar gas number density 



© RAS, MNRAS 000, 1-9 



6 Raffaella Schneider, Kazuyuki Omukai, Akio Inoue & Andrea Ferrara 



10^ 



10" 



10= 



102 



10 



: io"M,, 










~ PISN igSMg 










E [Z/H]= 


-~, -6, 


-4. -2 




: 


- [Z/H] = 


-7, -5, 


-3. -1 -~ 




- 


r p\ 










~ 1 W ■ y' 

\ ' \ 

\ \ \ 


\ 


/ \ . / 




- 


r :,A \^ 










\ 










■ \ 




















5 




10 15 


20 





log n„ (cm-3) 



10^ 



10" 



10=" 



10= 



10 



E ' 1 ' ' ' ' 1 ' ' ' ' 1 


1 1 1 1 1 1 


I 


- SN SSMg /' 




# 


E [7/H] = -». ' -fi -^ 




E 


[Z/H] = -7, -5, -3, 




■ - 


-fx ■ ^ 




— 


\ L\ .^^'^'^^^ 




\ 


" \ \ ^ ^^^^^ 


-5 // 






^ ^ y 
/ ,■ 




' \-2 






/ 








, , , 1 , , , 


1 = 



5 10 15 

log n„ (cm-3) 



20 



Figure 3. Thermal evolution of collapsing pre-stellar gas clouds as a function of the central gas number density enriched to different 
metallicities by the explosion of a PISN 195 Mq (left) and of a zero metallicity 22 SNII (right). Solid curves indicate total metallicities 
(gas-phase and solid dust grains) [Z/'H]=hog{Z /Zq) = — oo, —6, —4, —2 and dashed curves to [Z/H] = -7, -5, -3, -1. The diagonal dotted 
curves identify values of the Jeans mass corresponding to given thermal states. 
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Figure 4. Gas number density where 7 = 1 as a function of total metallicity (gas-phase metals and solid grains), for pre-stellar clouds 
enriched by the explosion of a 195 Mq PISN (left panel) and a 22 Mq SNII (right panel). Points indicated with circles (triangles) 
represent fragmentation induced by line (dust) cooling. Crosses show the effect of 3-body H2 formation and H2O cooling and, at the 
lowest metallicities, squares indicate fragmentation due to H2 collision-induced emission (see text). 



where this condition is met as a function of the initial total 
metallicity, for pre-stellar clouds enriched by the explosion 
of a 195Mq PISN (left panel) and zero metalhcity 22 Mq 
SNII (right panel). The "7 = 1" points can be classified 
into three main branches: the lowest-density points, shown 
as open circles, indicate fragmentation epochs which follow 
the temperature dip induced by line-cooling (hereafter line- 
induced fragmentation); the highest-density points represent 
fragmentation following cooling due to H2 collision-induced 
emission at the lowermost metallicities (squares) and due 
to dust thermal emission (hereafter dust-induced fragmen- 
tation) at higher metallicities (triangles); for completeness, 
we also show a third branch of 7 = 1 points (crosses) which 
corresponds to gas number densities rifrag ~ 10^" cm~^, and 
which is associated to cooling by 3-body H2 formation and 



H2O. As it can be inferred from Fig. 3, the cooling rate 
at these intermediate densities leads to temperature decre- 
ments which tend to be short-lived and shallow and are 
unlikely to induce fragmentation. The local maxima of the 
line-induced fragmentation branch result from the compe- 
tition between HD cooling (which starts to be effective at 
Z Si 10~^Zq) and heating due to H2 formation on grain sur- 
face which, for Z ^ 1Q~'^Zq prevents the gas from reaching 
the NLTE-LTE transition for HD at densities n ~ 10^ cm"^. 

It is important to stress that the presence of a temper- 
ature dip in the thermal evolutionary track as well as the 
condition that 7 < 1 imply only the possibility of fragmenta- 
tion. For clouds of solar metallicity, recent three-dimensional 
hydrodynamical simulations show that the single tempera- 
ture minimum present at a density comparable to that of 
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observed cloud cores plays an important role in determining 
the peak mass of the IMF (Jappsen et al. 2005). In contrast, 
the second minimum in the thermal evolutionary track due 
to dust-cooling occurs at a density which is several order 
of magnitudes higher, deep in the interior of the collapsing 
cloud. Therefore, it is not clear whether the same argument 
may be applied as well as the fraction of mass which may 
be affected. 

Dust-induced fragmentation will occur if very flattened 
or elongated structures form in the dense central region, 
and if the fragments formed do not subsequently merge. 
Indeed, following the first line-induced fragmentation, pre- 
stellar gas clouds experience runaway collapse, during which 
density perturbations might bo erased by pressure forces. If 
so, with little elongation of the densest core, dust-induced 
fragmentation at higher densities would be prohibited. By 
means of a semi-analytic model for the linear evolution of 
non-spherical deformations, Omukai et al. (2005) discuss the 
fragmentation histories of collapsing pre-stellar clouds and 
confirm that dust-induced fragmentation does occur even 
for metallicities as low as Z = 10"®^©. The mass frac- 
tion in low-mass fragments is initially small but becomes 
dominant for Z ^ W~^Zq and continues to grow with in- 
creasing Z. More recently, Tsuribe & Omukai (2006) used 
three-dimensional hydrodynamical simulations to follow the 
evolution of low-metallicity cores during the dust-cooling 
phase and discuss the conditions for fragmentation into low- 
mass clumps. Their results confirm that the previous anal- 
ysis of Omukai et al. (2005) was broadly correct: fragmen- 
tation does occur in the dust-cooling phase and sub-solar 
mass fragments are indeed produced. 

3.4 Characteristic mass 

In the present context, where we assume that thermal pres- 
sure is the main force opposing gravity, the characteristic 
mass of each fragment is given by the Jeans mass at frag- 
mentation epoch, 

Mfrag = Mj (nf rag, Tfrag) OC pX^j OC Tj^ag'^frag^- 

Fig. 5 shows the .Jeans masses associated to the lower 
and upper density branches represented in Fig. 4, hence 
to line-induced and dust-induced fragmentation. Since line- 
induced fragmentation occurs at relatively low densities, 
the corresponding fragment masses are large. In particu- 
lar, at the lowest metallicities we find that the character- 
istic fragmentation mass scale, induced by H2 line cooling, 
is Mfrag ~ 10^ M0, consistent with the results of numer- 
ical models. Still, even for total metallicities in the range 
10~^ < Z/Zq < 10~^, metal- line cooling leads to fragment 
masses M^ag > lOOM©. 

Conversely, dust-induced fragmentation leads to solar 
or sub-solar mass fragments, with O.OIM© < Mfrag < IMq. 
Indeed, this fragmentation mode occurs at much higher 
densities (see Fig. 4) where the associated Jeans masses 
are correspondingly lower. The characteristic fragment mass 
increases for increasing metallicity reflecting the opposite 
trend in the dust-induced fragmentation density shown in 
Fig. 4. For both supernova models, this fragmentation mode 
is active already at very low metallicities and the solid line 
illustrates the transition in fragmentation mass scales for 

Zcr = lO-^Z©. 



The importance of dust for the formation of low-mass 
stars was first emphasized by the pioneering work of Low 
&i Lynden-Bell (1976), where they also recognized that dust 
mantains its cooling efficiency down to a reduction by mass 
of 10"'"' of the present-day value. 

The results that we find suggest that in the absence 
of dust pre-stellar clouds enriched by primordial supernovae 
are not able to form solar or sub-solar mass fragments. 



4 DISCUSSION 

In this paper we have presented the first self-consistent 
model of thermal evolution of pre-stellar gas clouds en- 
riched with metals and dust by primordial supernovae. We 
have focused on two representative cases: a 195 M© pair- 
instability supernova and a 22 M© Typell supernova. Us- 
ing the one-zone semi-analytic model initially developed by 
Omukai (2000), and recently improved in Omukai et al. 
(2005), we have followed the thermal and fragmentation his- 
tory of pre-stellar clouds enriched by the products of single 
primordial supernovae. We have sampled initial metallicities 
in the range < Z/Zq < 0.1, as representative of different 
enrichment histories, which may vary depending on metal 
diffusion and mixing. 

Neglecting differential enrichment by metals and dust 
grains, the initial composition of the collapsing pre-stellar 
gas clouds directly reflects the nucleosynthetic products of 
the polluting SN, its specific dust depletion factor, rela- 
tive gas-phase metal abundances and dust grain species and 
sizes. We have therefore modified the model to allow a self- 
consistent description of metal-line cooling, dust thermal 
emission, dust-gas collision rate, H2 formation rate on dust 
grains. In all theoretical studies carried out so far, these 
processes were implemented under the implicit assumption 
that the gas-phase metal abundances and dust grain proper- 
ties were the same as those observed in the local interstellar 
medium. The only exception is the study of Bromm & Loeb 
(2003), where separate critical abundances of CII and 01 
enabling cooling and fragmentation in collapsing gas clouds 
were derived. Still, they limited their study to the thermal 
evolution in the low-density regime, including only metal 
line cooling and fragmentation. 

The main results of our study can be summarized as 
follows: 

(i) As long as thermal physics is considered to play a ma- 
jor role in setting a characteristic stellar mass scale, dust- 
induced fragmentation appears to be the most promising 
mechanism for the formation of solar and sub-solar mass 
fragments. Since dust can be promptly synthesized in the 
ejecta of primordial SN, it is likely that its contribution to 
thermal and fragmentation history of collapsing clouds at 
high redshifts had been relevant. 

(ii) Due to the strong depiction of C atoms onto AC 
grains and CO molecules for most SN models, our analy- 
sis shows that CII makes a negligible contribution to fine- 
structure line cooling, leaving OI as the main cooling agent 
in the gas phase. 

(iii) Line-cooling, being limited to the initial evolu- 
tion of collapsing gas clouds, for densities in the range 
10* - 10'^ cm , induces fragmentation into relatively 
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Figure 5. Thermal Jeans mass corresponding to line-induced fragmentation (circles in Fig. 4) and dust-induced fragmentation (triangles 
in Fig. 4) as a function of total metallicity, for pre-stellar clouds enriched by the explosion of a 195 Mq PISN (left panel) and a 22 Mq 
SNII (right panel). The solid line indicates the transition in characteristic fragment masses when the dust-induced fragmentation mode 
is activated at a critical metallicity Zct = 1Q~^Zq. 



large clumps, with characteristic fragment masses Mfrag > 
IOOM0, even for initial gas metalhcitics Z ^ IQ^^Zq. 

(iv) The formation of solar or sub-solar mass fragments 
might operate already at metanicities Zor = 1Q~^Zq. The 
results that we find suggest a re-formulation of the critical 
metallicity scenario: indeed, once a self-consistent model is 
adopted, the presence of gas-phase metals can not by it- 
self lead to solar or sub-solar characteristic masses whereas 
primordial supernovae of all masses (that is both in the stan- 
dard SNII mass range and in massive PISN), appear to de- 
plete enough metals in solid grains to enable the formation 
of low-mass fragments. 

Our approach to determine the minimum fragment 
masses as the Jeans masses corresponding to the last frag- 
mentation (temperature dip and 7 = 1 condition) may 
be an oversimplification of the complex dynamical evolu- 
tion of cloud cores. Still, our findings arc supported by re- 
sults of a semi-analytical model for the linear evolution of 
non-spherical deformations by Omukai et al. (2005) and by 
three-dimensional hydrodynamical simulations of Tsuribe & 
Omukai (2006), which confirm that fragmentation occurs 
during the dust-cooling phase and that sub-solar mass frag- 
ments are indeed produced. 

Let us emphasize that other processes might be equally 
important in reducing the characteristic mass scales of 
metal-free or very metal-poor star forming regions. Among 
these, the bimodal IMF suggested by Umemura & Nakamura 
(2002), the presence of a strong UV field in the vicinity of 
a very massive star (Omukai & Yoshii 2003), HD cooling 
triggered in relic HII regions around very massive stars af- 
ter the death of the exciting star (Uehara & Inutsuka 2000; 
Nagakura & Omukai 2005; Johnson & Bromm 2006), shock- 
compression induced by the first SN explosions (MacKey, 
Bromm & Hernquist 2003; Salvaterra, Ferrara & Schneider 
2004). 

Observationally, stellar relics in the halo of our Galaxy 
represent fundamental records of low-mass star formation at 
high redshifts. The HK-survey, Hamburg-ESO Prism Sur- 



vey and Sloan Digital Sky Survey provide already a sample 
of about 5000 stars with [Fe/H]< —2, among which 750 
with [Fe/H]< -3. The project SEGUE (Sloan Extention 
for Galaxy Understanding and Exploration) shoud identify 
about 25000 stars with [Fe/H]< -2 (Beers & Christlieb 
2005). The interpretation of shape and low- metallicity tail 
of the metallicity distribution function shoud provide impor- 
tant insights on the origin of the first low-mass stars that 
have formed in the Universe (Schneider et al. 2003; Tumlin- 
son 2006; Salvadori, Schneider & Ferrara 2006). 

The surface elemental abundances of the observed stel- 
lar relics should shed light on the IMF of the stars mostly 
responsible for early metal enrichment (Tumlinson, Venkate- 
san & Shull 2004; Iwamoto ct al. 2005; Beers & Christlieb 
2005). In principle, dust depletion docs not affect the inter- 
pretation of the observed chemical abundances because, once 
sufiiciently high central densities axe reached, dust grains 
present in collapsing pre-stellar clouds are destroyed and 
the original gas-phase metal abundances (reflecting the nu- 
cleosynthetic products of the SN progenitor star) are re- 
established. An alternative scenario has been recently pro- 
posed by Venkatesan, Nath & Shull (2006), who suggest that 
the elements forming the composition of the dominant dust 
compounds created in primordial supernovae (Si, Mg, O, 
C) can be selectively transported in the UV radiation field 
from the first stars and decoupled from the background SN 
metals in the gas phase. Interestingly, these elements appear 
to have enhanced abundances in the most metal poor halo 
stars observed. 
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